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SECURITY  I 


A  MECHANISM  FOR  HEATING  THE  SOLAR  CORONA 

The  nature  of  the  mechanism  which  causes  temperatures  in  the 
solar  atmosphere  to  rise  from  about  5000  K  in  the  photosphere  to 
over  10^  K  in  the  corona  is  a  major  unsolved  problem.  Wave  processes 
have  been  invoked  to  transport  the  necessary  energy  from  the  sun's 
surface  (Osterbrock,  1961).  Ohmic  heating  or  magnetic  reconnection 
has  also  been  proposed  (Vaiana  and  Rosner,  1978).  It  is  the  purpose 
of  this  letter  to  point  out  that  another  process,  compress ional  work 
done  on  the  plasma  by  expanding  magnetic  flux,  may  be  important. 

Observations  made  in  the  last  few  years  (Sheeley,  et  al.,  1975) 
show  that  the  corona  is  highly  structured,  with  magnetic  fields  playing 
a  signf leant  -  perhaps  dominant  -  role.  They  also  reveal  that  the  corona  is 
in  a  continual  state  of  fluctuation  on  both  large  and  small  space-time 
scales  (Wentzel,  1978).  In  spite  of  this,  elevated  coronal  temperatures 
are  apparently  maintained  continuously.  It  is  thus  possible  to  model 
the  coronal  heating  as  the  outcome  of  frequently  occurring  but  discrete 
events,  each  of  which  deposits  some  fraction  of  the  energy  needed  to 
balance  the  power  dissipated  by  the  corona. 

Consider  a  bubble-like  object,  composed  of  magnetic  flux  and  plasma, 

emerging  from  the  surface  of  the  sun.  The  exact  topology  of  this 

entity  is  not  critical;  it  may  resemble  a  jet  more  than  a  bubble,  for 

example,  or  it  may  have  the  form  of  an  arch,  and  the  field  lines  may 

remain  connected  to  the  solar  magnetic  dipole.  The  plasma  density  is 

less  than  or  equal  to  that  in  the  surrounding  medium,  but  the  Initial 

2 

magnetic  pressure  B  /8k  Is  assumed  to  exceed  the  average  ambient  pressure, 

so  that  the  structure  expands  as  it  rises.  This  continues  until  pressure 

balance  is  reached,  reconnection  annihilates  the  magnetic  flux,  or  the 

bubble  vents  out  the  top  of  the  atmosphere.  In  any  case  the  work  AW 

2 

done  on  the  surrounding  plasma  is  essentially  AVB  /8k,  where  AV  is 
the  initial  volume  of  the  bubble.  Moreover,  a  mass  Am  •  pAV  is  injected 
into  the  corona  at  the  same  time,  where  P  is  the  masa  density  when  the 
bubble  emerges.  (  \  j  |  1  1  I 

l  ml 
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If  structures  of  this  type  are  continuously  generated  so  that 
the  average  volume  of  flux  emitted  per  unit  time  is  V  and  if  they 
are  responsible  for  heating  the  corona,  then  power  balance  implies 
that  in  the  quiet  sun  (Withbroe  and  Noyes,  1977)  ,  on  the  average. 


V  B2/8ir 


■/i p  ,  +  p  +  p 

l  rad  conv  condl 


t*  (105  +  2  x  105  +  5  x  104)  A  erg  s-1 


■  3.5  x  105  erg  s  1.  (1) 

2 

Here  A#  ■  4ttR#  and  the  three  terms  on  the  right  hand  side  of  (1)  represent 
losses  due  to  radiation,  convection  (i.e.,  the  solar  wind),  and  thermal 
conduction  along  the  field  lines  to  the  cooler  chromosphere.  The 
structures  emanate  with  a  velocity  v  from  some  fraction  F  of  the  sun's  bi- 
polar  magnetic  regions.  The  latter  have  an  area~0.1A  ,  so  V~0.1  F  A  v. 

If  the  emerging  structure  has  low  density  and  rises  because  of  buoyancy 
forces,  then  v  ~  (gtQ)  ,  where  g  is  the  surface  gravitational  acceleration 
and  Tq  is  a  characteristic  dimension  for  bipolar  magnetic  regions.  If 
Instead  it  is  propelled  by  magnetic  forces,  then  v  approximately  equals 
the  Alfven  speed  »  B  (4np)  .  Other  mechanisms  lead  to  different 

estimates  of  v.  Although  the  choice  is  not  crucial  to  the  theory,  we 

will  assume  here  that  v  »  v,  ~  2  x  10^  B  cm  s  1  (with  B  in  G),  where 

A  10  -3 

p  has  been  calculated  by  taking  an  ion  number  density  ~  10  cm  .  Then 

from  (1)  a  fraction  F~40  B  of  the  sun's  magnetically  active  regions 

must  be  emitting  flux  at  any  Instant  of  time.  The  requirement  F  £  1  implies 

that  B  £  3  G  must  hold. 

Similarly,  the  upward  motion  produces  a  surge  which  injects  mass  Into 
the  corona  at  a  rate 


*  *  -7  -2  -1 

a  -  p  V  ~  1.3  X  10  B  g  s  . 


2- 


The  requirement  that  the  velocity  of  the  material  sprayed  out  of  the 

solar  atmosphere  not  exceed  the  sum  of  the  observed  solar  wind  speed 

and  the  escape  speed  implies  v  &  10  cm  a,  or  B  <  50  G.  At  this  upper 

-11  -1 

limit,  the  total  mass  injection  rate  is  m  ~  4  x  10  A(  g  8  .  This 

is  comparable  with  observed  solar  wind  mass  losses  (Withbroe  and  Noyes, 
1977) . 

To  obtain  more  information  about  how  an  expanding  magnetic  "bubble" 
heats  the  corona  as  a  function  of,  e.g.  altitude,  we  need  to  adopt  a 
specific  model.  The  following  discussion  illustrates  this.  The  model 
chosen,  although  unrealistic  in  some  respects,  makes  it  clear  that  wave 
and  current  dissipation  are  not  required  by  the  process.  Consider  the 
time-dependent  axisymmetric  system  given  by 


po f"9; 

(3) 

a!  p  «2/t2; 

(4) 

r/t; 

(5) 

z/4t; 

(6) 

-  4>  (r/L)  sin 

(2z/L)/u02t2; 

(7) 

■  $  (r/L)  cos 

(2z/L)/u02t2; 

(8) 

■  $  cos  <2z/L)/u02t2. 

(9) 

Here  r  and  z  are  distances  measured  from  and  along  the  axis  of  symmetry; 
u 

f  *  (t/tg)  ;  L  *  Xf;  and  pg,  Ug,  tg,  1  and  $  are  constants.  Equations  (3)- 
(9),  which  can  readily  be  shown  to  satisfy  the  ideal  MHD  equations  with 
adiabatic  index  y  *  5/3,  describe  a  force-free  self-similar  configuration 
with  similarity  variables  £  ■  r/ugt  and  c  ■  2z/L,  corresponding  to  independent 
uniform  homologous  expansions  in  the  transverse  and  vertical  directions 
with  velocities  u  and  v,  respectively.  All  three  components  of  the  magnetic 
Induction  £  are  proportional  to  the  flux  *  contained  within  a  tube  of  radius 
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u 

R(z,t)  "  i^t  sec  (2z/L),  and  hence  can  be  chosen  arbitrarily  large 

without  affecting  the  material  praaaura  p,  the  density  or  flow  field. 

Equations  (3)  -  (9)  become  singular  at  t  ■  o,  and  so  must  be 

initialised  at  some  finite  time.  If  we  take  tQ  to  be  the  initial  time, 

then  p0,  rQ  •  uQt0,  A,  and  1Q  »  6/r02  can  be  Interpreted  as  the  initial 

values  of  the  density,  characteristic  horlsontal  and  vertical  length  scales, 

and  characteristic  magnetic  field  strength.  Likewise,  the  flux  tube  radius, 

and  therefore  and  on  the  flux  tube  boundary,  diverge  at  v.  -  *L/4, 

so  the  solution  must  be  restricted  to  a  region  0  <_  s  <  Z  ■  zQf,  where 

zQ  is  a  constant  (to  be  specified  below)  satisfying  zQ  <  *d/4.  For 

z  >  Z  the  system  is  not  descrlbable  by  a  similarity  solution,  or  it  must 

violate  the  assumption  of  ideal  force-free  flow. 

Ue  will  overlook  the  problems  with  the  magnetic  field  above  z  ■  Z, 

as  well  as  the  neglect  In  the  model  of  the  sun's  curvature  and  gravity, 

and  will  Introduce  the  further  simplifying  assumption  that  the  plasma  in  the 

exterior  region  r  >  R  is  unmagnetized.  At  the  flux  tube  boundary  a  kinematic 

2  2 

condition  and  the  condition  of  pressure  balance,  p#jjt  ■  p  +  B  /gw  s»  B  /8w, 

must  be  satisfied,  where  in  consequence  of  our  motivating  physical  picture 

2 

we  ignore  p  compared  with  B  /8v.  It  is  unnecessary  to  construct  a  complete 

model  of  the  flow  in  the  exterior  region.  Instead,  we  limit  ourselves  to 

calculating  the  work  done  by  the  expanding  flux  tube  on  the  surrounding 
medium. 

The  power  per  unit  height  generated  by  the  expansion  Is 

2  ig 

“  2r>  R  1 1  +  (||)  1  ft  •  X  P.xt  "  J  I  u(R)-v(R)  B2 

•  (t2 /Ug2t^)  (  1  -  c/8  tan  (.  }  f  (Ugt/L)2  ♦  cos^c]  sec2  t,  (10) 

where  v  is  the  velocity  of  the  flux  tube  boundary.  He  note  that  (10) 

becomes  negative  for  C  sufficiently  cloae  to  «/2 .  There  the  tube,  which 

flares  out  almost  horizontally,  ia  locally  receding  from  the  plasma.  To 

avoid  this,  we  choose  z  so  that  (z  /A)  tan  (2s  /A)  -  4,  i.e.,  s  m  0.699*. 

0  0  0  0 
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Me  will  mw»  that,  whatever  tha  actual  fora  of  eha  configuration  for 
•  >  Z,  tha  exterior  medium  doe*  aa  Much  work  on  It  aa  It  does  on  tha 
medium.  Than  If  tha  expansion  work  dona  at  a  given  altitude  roawina 
thara  in  tha  fora  of  tharaal  energy,  tha  haat  deposited  aa  a  function  of  z  is 

glvan  by  «  (a)  •/  "  dt  •  In  ri*‘  l»  #  *  W  a  r0  Plottad  ogoinnt 

*0  °* 

a  for  thraa  different  values  of  tha  aspect  ratio  a  ■  d/rft.  It  can  be 

2  2-2  0 

seen  that  H  scales  roughly  aa  r^  Eg  a  For  a  ~  1  it  has  a  maximum  at 

a  ~  M  ••  ®  •*  •  •  tha  heating  becomes  strongest  at  tha  bo t ton  of  tha 

structure.  In  either  case,  nost  of  the  haat  is  depoaltad  in  a  distance 

of  order  rQ  (tha  initial  horizontal  scale  of  the  bipolar  magnetic  region*. ) 

fro*  tha  bottoa  of  tha  atmosphere.  Integration  of  (10)  over  both  a  and  t 

shorn  that  the  total  heat  evolved  in  the  expansion  process  la  equal  to 

the  magnetic  energy  Initially  contained  in  the  configuration  (pro|M*ctlonal 
2  2 

to  Eg  rg  X),  aa  it  should.  The  actual  manner  in  which  the  unmagnetized 
atmosphere  is  heated  depends  on  convection,  conduction  and  radiation  pro¬ 
cesses,  and  Is  of  course  not  local.  Evidently  tha  rapid  expansion  of 
the  magnetic  bubble  can  drive  shocks,  particularly  if  the  expansion 
velocities  are  as  large  as  v^.  Numerical  models  are  currently  being 
developed  at  the  Naval  Research  Laboratory  in  order  to  make  global 
quantitative  predictions  of  the  heating  to  be  expected  from  the  mechanism 
discussed  here,  together  with  the  effects  of  wave  and  ohmic  heating  and 
other  plasma  processes. 

In  suasHry,  bubbles  or  other  structures  containing  magnetic  flux 

with  field  strengths  ~  20C  and  having  transverse  dimensions  like  those 

of  the  initial  stages  of  bipolar  magnetic  regions,  can  supply  heat  to  the 

solar  atmosphere  through  expansion  in  amounts  sufficient  to  explain  the 

observed  elevated  coronal  temperatures.  If  magnetic  forces  are  responsible 

for  expelling  this  flux,  the  bubbles  should  be  emerging  from  a  fraction 

F  >  0.03Z  of  the  bipolar  magnetic  regions  at  any  time,  with  characteristic 
“  _1  3 

speeds  v.  <  500  km  s  and  lifetimes  R  /v.  >  10  a,  and  should  deposit 

A  +*  o  A  45 

their  energy  principally  within  vertical  distances  ~  10  -  10  km  from 

the  surface  of  the  sun. 
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Fig.  1  —  Reduced  heating  function  R  vs.  altitude  for  three  values  of  the  aspect  parameter 
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